Thermodynamics of the half-filled Kondo lattice model around the atomic limit 
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We present a perturbation theory for studying thermodynamic properties of the Kondo spin liquid 
phase of the half-filled Kondo lattice model. The grand partition function is derived to calculate 
chemical potential, spin and charge susceptibilities and specific heat. The treatment is applicable to 
the model with strong couplings in any dimensions (one, two and three dimensions). The chemical 
potential equals zero at any temperatures, satisfying the requirement of the particle-hole symmetry. 
Thermally activated behaviors of the spin(charge) susceptibility due to the spin(quasiparticle) gap 
can be seen and the two-peak structure of the specific heat is obtained. The same treatment to the 
periodic Anderson model around atomic limit is also briefiy discussed. 
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PACS numbers: 75.10.Jm, 71.27.+a, 71.10.Fd 

The Kondo lattice model (KLM) and the periodic An- 
derson model (PAM) are two prototype models for the 
heavy fermion system and related materials |[ ^ . Both 
models describe a band of itinerant conduction electrons 
interacting with a lattice of magnetic impurities. At half- 
filling, they are usually viewed as standard models to 
understand the so-called Kondo insulator. The PAM is 
written as 

{i,j),s i,s 



(1) 



where the fermion operators cl^ if is) create conduction 
(impurity) electrons and the summation j) taken 
over all nearest neighbors. In the Kondo limit, Uf » 
V and ep ^ £/ (ef is the Fermi energy), the charge 
fluctuation of / electrons is suppressed and each impurity 
site is occupied by one and only one / electron. Then the 
PAM can be simplified and reduced to the KLMB, 



KLM 



H.c. 



(2) 



with the spin-exchange coupling J = 8V'^/Uf under 

symmetric conditions. Here s'^ — J2s s' l'^s,s'clgCis' is 
the spin-density operator of conduction electrons and 
S = J2s s' \'^s,s' fisfis' denotes the localized spin with 
fis if is) being the pseudo- fermion operators. (Ts^s' is the 
Pauli matrix. 

The purpose of this work is mainly to consider the half- 
filled KLM whose ground state and finite-temperature 
properties have been intensively studied in recent years [|[ 
1' III 1' 0' "^^^ ground-state phase diagram has been 
now well established Owing to the Kondo screen- 
ing effect, the one dimensional (ID) model has a Kondo 
spin liquid (KSL) state at all values of the coupling J/t. 
In higher dimensions, it is suggested that competition 
between the Kondo screening and the Ruderman-Kittel- 
Kasuya-Yosida interaction leads to quantum phase tran- 
sitions from the KSL state to the antiferromagnetically 



ordered state as J/t decreases, with the critical point 
{J/t)c ~ 1.4 for two-dimensional (2D) and 2.0 for three- 
dimensional (3D) KLM'sg 0, |]. The KSL phase is 
characterized by three different energy scales: the spin 
(As), charge (Ac) and quasiparticle (A^p) gaps. The 
spin gap vanishes at the critical point while the other 
two gaps still remain finite in the ordered phase. The 
finite-temperature properties of the KSL phase are dom- 
inated by all three energy scales. A number of power- 
ful numerical methods including quantum Monte Carlo 
(QMC) simulationsB] , the density matrix renormaliza- 
tion group method ff and the finite-temperature Lanczos 
technique |lo[] have been used to investigate the thermo- 
dynamic properties of the KLM in ID and 2D. However, 
less analytical work is known, especially for the higher 
dimensional models. 

In this paper, we propose a finite-temperature pertur- 
bation theory to study the thermodynamics of the KSL 
phase around the atomic limit t ^ 0. We believe that 
the approach is applicable to the KLM in any dimensions 
and at any temperatures as long as the coupling is suf- 
ficiently strong, J ^ t. The model under consideration 
is a generalized KLM by introducing an on-site Coulomb 
interaction between conduction electrons to the Hamil- 
tonian (2), 



H = H 



KLM 



(3) 



The hopping term in Eq. (2) is treated as a pertur- 
bation. A similar treatment has been applied to pure 
spin systems with the dimerized ground state and 
yields good agreement with experiment results and QMC 
simulations . 

The strong coupling limit has been taken as a good 
starting point to understand the ground state properties 
of the KSL phase. Various techniques based on this limit, 
such as direct strong-coupling expansions!^, slave par- 
ticle approach fl^, [l^ and projector techniques [p^, are 
employed to calculate the ground state energy and the 
excitation spectra, and even to determine the quantum 
transition point of higher-dimensional KLM's. In the 
limit case, the ground state and low-lying excited states 
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TABLE I: Eigen states and energies of each atom, h is the 
external magnetic field. = 1,0, —1 is the spin-z component 
of triplet excitations and s = | ( — |) corresponds to spin up 
(spin down) of hole and double occupancy states. 
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can be described by simple wave functions and thus it 
is much easier to define the energy gapsQ. For each 
atom, the Hilbert space consists of eight possible quan- 
tum states (see table I): a singlet state |si} = — 

cl^f-^)\Q), three degenerate triplet states |4) — '^Tfi'jl^)' 

\t7') = 4ifl\o) and = ^{cyi + cyi^m, two 

hole states \ais) — fjjO) and two double occupancy states 



\bis) = C'T^i/islO)- The unperturbed lattice eigen states 
are simply the direct products of the atomic states, e.g., 
\siS2---Sn) is the ground state. The lowest spin excitation 
is promoting one singlet to triplet. And the lowest charge 
excitation is to create a pair of one vacant and one dou- 
bly occupied site, with a spin-singlet combination of the 
two localized spins. The quasiparticle is defined as either 
one vacancy or one double-occupancy configuration. One 
can get a qualitatively correct picture of the KSL phase 
for smaller couplings from the limit case. Naturally, we 
expect that it is also an appropriate starting point for 
understanding thermodynamic properties. 

We use the grand canonical ensemble. First the 
grand partition function is calculated, Z((3,fi) = 
Tr[exp{—f3{H — fin''))] where f3 = l/fc^T, is the chemi- 



J 



cal potential and n'^ is the number operator of conduction 
electrons. Since the / electrons are completely localized 
and they have no contributions to the fluctuation of the 
electron number in the KLM, only the number of con- 
duction electrons is considered. Performing the cumulant 
expansion, the partition function reads 



oo 



n=l 



(4) 



with the cluster function C/„ defined as 



^0 



Zq 2 \Zo 



Us 
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Here Z„ is the n-th order perturbation expansion of the 
partition function. 



Z„ = (-1)" y"dTl J dT2... 



X J dT,,Tr[e'^^"°->^^^^H\T^)H\T2)...H\T,,)l<6) 



and H\t) = e-r(H°-i,n-)jjI^-r{H°-i,n-)^ jjg^.^ 

be integrated out in the basic set of eigenstates of the 
unperturbed Hamiltonian |pT| . It contains the product of 
n perturbation matrices. 
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{Hl^)mn = -t{m\ I]s(4,sCj,s + H.c.)\n) is the matrix el- 
ement of the perturbation in the Hilbert space of the 
atom lattice. Nevertheless only those terms in which all 
the sites ii,ji;«2,i2; ■■■;*n, jn are connected by the inter- 
action are taken into account in the cluster function Un- 
The computation of the perturbation expansions is 
straightforward but very laborious. Above all, one has 
to compute the perturbation matrix elements {Hlj)mn- 
It is important to note that both fg and Cg are fermionic 
operators and the sign changes when they exchange their 
sites. For example. 



(7) 



\ais) and \bis) are fermionlike states, while \tf) and \si) 
are bosonlike states, 



\0:isbjs' ) 



-\bjs'ais), 



(8) 



{Hlj)mn is a real, diagonal matrix. For each nonzero 
element, the final and initial states (m| and \n) differ 
only in the states of the two adjacent atoms and the rest 
of the lattice is unchanged. So only 8x8 different states 
are concerned. The interaction term just transfers an 
electron from one site to its neighbors and so it makes 
the states of the two involved sites definitely changed but 
conserves the number of electrons and the z component 
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FIG. 1: Spin (a) and charge (b) susceptibilities Xs, Xc vs T jt 
for the ID KLM with = 0. 
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FIG. 2: Specific heat C{T) vs T/t of the ID KLM with 
(a) (7 = and J/t = 4.0,3.5,3.0,2.5; (b) J/t = 2.5 and 
Uc = 0,1.0,2.5. 



of spins. Taking advantage of this point, the calculation 
of the matrix elements and the cluster functions can be 
greatly simplified. For example, one can easily get that 
all diagonal elements of perturbation matrix are zero. 
To the d-dimensional cubic lattice KLM, all the odd- 
order cluster functions also vanish. In this paper, we 
only evaluate the second-order expansion [n = 2). 

Before calculating thermodynamic properties, one 
must first determine the chemical potential ji by solving 
the equation 



(9) 



where (n'^) is the average number of conduction elec- 
tions per site, N is the number of sites and F{(3,^) — 
— ■g/n[Z(/?, /i)] is the free energy. Generally, the chemical 
potential is a function of the temperature. Nevertheless, 
/i should be exactly equal to zero at any temperatures 
due to the particle-hole symmetry at half filling. There 
is one conduction electron per site on average in this case. 
Setting {n'^) = 1 in the above equation, we obtain that 
the requirement is well satisfied up to the second-order 
expansion. 

We proceed to calculate the spin susceptibility Xs [T) — 
-j^{%f)h=o, charge susceptibility Xc{T) = and 

specific heat C{T) = -jfP'^^is^- AU the thermody- 
namic quantities are derived from the free energy by a 
standard procedure. Fig. (1) shows numerical results of 
spin and charge susceptibilities for the ID KLM. Raising 
temperature from zero, both XsiT) and XciT) rise expo- 
nentially, indicating the existence of the spin and quasi- 
particle gaps, and then they reach a maximum. The 



peaks of Xs{T) and XciT) move to lower temperatures 
with J/t decreasing, suggesting that the A.; and Aqp di- 
minish correspondingly. As discussed by Shibata et al. |^ 
and Haule et al||l^, both As and Aqp are responsible 
for the spin susceptibility but mainly the lower one dom- 
inates the low temperature behaviors. In contrast, Xc is 
governed by Aqp only. Since Aqp < As near the atomic 
limit, both Xs{T) and Xc{T) are dominated by the quasi- 
particle gap. We find that the activation energies esti- 
mated by fitting the spin and charge susceptibilities are 
similar when J/t > 4, consistent with the conclusion of 
the numerical results l^], |l^. However the second-order 
perturbation result is not accurate enough to discuss the 
delicate relation between the activation energies of Xs (T) 
and Xc(T) and the spin and quasiparticle gaps for smaller 
couplings. 

For J/t ^2.5, there appears an unphysical protrusion 
on the peak of Xc (T) . This value is an approximate lower 
limit of the J/t range for which the perturbation theory is 
reasonable for all temperatures. Within the second-order 
perturbation theory, the lower limit of J/t is about 2.5, 
3.0 and 4.0 for ID, 2D and 3D KLM's respectively. This 
is consistent with the zero-temperature strong-coupling 
expansions by which the spin gap within second order is 

(2) 

As ' tends to zero quickly at 



A (2) 
given as As 



J- 



20<ir 
3J 



J/t w 2.58, 3.65 and 4.47 for d = 1, 2 and 3. Considering 
higher-order expansions is expected not only to improve 
the accuracy of the result, also to reduce the lower limit 
of J/t. This task might be achieved by the high order 
series expansion technique [l^ which has been successful 
in explaining ground-state properties of KLM[Q. The 
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FIG. 3: Specific heat of the 3D KLM. 



methodology for carrying out high order expansions at 
finite temperature for spin systems has been devefoped 
recently Q. 

Results of the specific heat for ID and 3D KLM's are 
presented in Figs. (2) and (3). The specific heat contains 
information on both spin and charge degrees of freedom. 
The charge fluctuation originates from the movement of 
conduction electrons. Moving an electron from one site 
to its neighbors on the ground state background creates 
a pair of one vacant and one doubly occupied site. The 
energy cost just is the charge gap. The charge fluctu- 
ation is greatly suppressed due to the large charge gap 
near the atomic limit, so only one peak attributed to the 
spin excitations can be seen. The location of the peak 
depends on the spin gap which varies with J. The hoping 
term stirs up the charge fluctuations. Then a new peak 
originating from the charge degree of freedom begins to 
be visible as J/t is reduced to 3.0, as shown in Fig. 2(a). 
The new peak corresponds to the specific heat of the free 
conduction electrons and is independent of J. The on-site 
Coulomb repulsion in the conduction electron band will 
enhance the charge gap and so it suppresses the charge 
fluctuation. As shown in Fig. 2(b), the second peak in 
C(T) will vanish with increasing Uc- Comparing Figs. 
(3) and (2), the results for the 3D KLM exhibit similar 



features to those for the ID case. As the T = per- 
turbation theory shows, various properties, including the 
characteristic energy gaps, of the KSL phase are maiiily 
determined by the local dimers in the atomic limit 
so the dimensionality does not make much difference in 
thermodynamic properties. This point is also indicated 
in numerical results for the 2D KLMjlOt. 

At last, we notice that the above treatment can be also 
applied to deal with the PAM, but with more complexity. 
Contrary to the KLM, the charge degrees of freedom of 
/ electrons must be considered in the PAM. The vacancy 
and double occupancy of the impurity site are possible 
and thus the Hilbert space of each atom is greatly en- 
larged. In general case each atom has 16 quantum states. 
The atomic limit of the PAM has been investigated by 
many authors |2^, |2l|. Choosing the atom Hamil- 
tonian as the unperturbed Hamiltonian, we can evaluate 
the thermodynamics near the atomic limit by perturba- 
tion theory as above. The physical properties of the PAM 
are more complicated than those of the KLM. For exam- 
ple, the specific heat will exhibit two peaks even in the 
atomic limit More recently, Moskalenko et ai. cal- 
culated the one-particle Green's function of the PAM by 
the hopping perturbation treatment [p3| . 

In summary, we have studied thermodynamic prop- 
erties of the Kondo spin liquid phase of the half-filled 
Kondo lattice model by a finite-temperature perturbation 
theory. The chemical potential, spin and charge suscep- 
tibilities and specific heat are calculated. The results are 
consistent qualitatively with those obtained by numeri- 
cal methods. It proves that the strong coupling limit is 
a reasonable starting point to study thermodynamics of 
the Kondo spin liquid phase. The accuracy of obtained 
results is expected to be improved by taking higher-order 
expansions into account. 

The author would like to thank P. Thalmeier for help- 
ful discussions. 



APPENDIX A 

The free energy per site within second-order perturba- 
tion is given as 
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4. 



where = 1, 2, 3 is the dimensionality and 



J 



zo = [l-Fe^^(^=-2^"'')l /(/i)-He-'^(-T-A'"'=) +e-'^(*-^"') [l + /(2/i)], 
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with /(/i) = e-^f^^' + eif^^. 
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